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abstract 



An effective algebraic approach to ^-matrix factorization into Jost 
matrices is developed in the case of coupled channels. The Jost matrix 
is given as a solution of boundary value Riemann - Hilbert problem. 
A rational form is assumed for tangent of mixing angle, while there is 
no limitations for approximation of phase shifts. 
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Introduction 



As it is well known, the Jost function plays an important role in scat- 
tering theory [jj], and particularly in the inverse scattering problem 
[)2|-[)5|. One of the methods of obtaining Jost function F{k) is to solve 
Riemann - Hilbert problem for the half-plane. In the frame of scat- 
tering theory the solution of this problem is reduced to the (S-matrix 
factorization in terms of Jost matrices The Jost matrix is also 
used to obtain the solutions of nonlinear differential equations (see, 
e.g. i, i). 

At present time a number of different approaches to the Jost ma- 
trix construction exists, e.g. [^-[[Mj. However, while the boundary 
value problem solution for one channel elastic scattering can always 
be represented explicitly in terms of phase shift, this is not so for 
two-channel scattering. In general case the Jost matrix can not be 
constructed effectively. It is known only that the matrix boundary 
value problem can be reduced to a system of singular integral equa- 
tions [T^ . While solving it one needs to take into account the following 
important fact. Jost matrix arises at an intermediate stage in the in- 
vestigation of various problems, so its explicit form must be the most 
adequate for analytical and numerical calculations. 

In this connection, one is interested to obtain the solution avoiding 
integral equations. 

Let us note that the procedure of obtaining the Jost matrix can 
be effectively simplified if the ^-matrix is factorized previously 0. 
This enables one to write the result in a compact form, as can be seen 
below from the concrete example. 

In this paper we shall use the approximate form for mixing param- 
eter. Some of preliminary results were presented in |TB1 . 



This paper is organized as follows. In Sec. || the boundary value 
Riemann - Hilbert problem for ^-matrix is formulated. The special 
features of the problem in matrix case are discussed in Sec. 0. Our 
algebraic approach consists of three steps. In Sec. fTI| the first step 
5'-matrix factorization — prefactorization is realized. The "correct- 
ing matrix" concept needed to solve the boundary value problem is 
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discussed in Sec. |V]. In Sec. |V] the final result for Jost matrix is ob- 
tained. The effectiveness of our approach is demonstrated in Sec. |V| 
for the well-known effective radius approximation. 



I The formulation of the problem 

The general mathematical formulation of the problem is the following. 
One needs to find the piecewise holomorphic matrix F{k) (i.e. its 
matrix elements are piecewise holomorphic functions) in the upper 
and lower half-planes. Boundary values of F{k) satisfy the following 
condition on the real axis: 

F+{k) = S'\k) F4k) , ImA: = 0, -oo < A: < oo (1) 

Here F^{k) is a matrix holomorphic in the upper half-plane, F^{k) — 
in the lower one, k is momentum in the center-of-mass system, S{k) 
is the scattering matrix which is non-singular: 

detS{k)y^O, lmk = (2) 

and Holder condition is valid for its matrix elements Sij{k) , i,j = 1,2 
on the real axis 

\s^J{kl) - s,jik2)\ < A\ki - /csl^ ^ > , < ^ < 1 , J = 1, 2 . (3) 

Standard conditions are imposed on the scattering matrix S{k) and 
on the F{k) matrix (the boundary values of the matrix F{k) are the 
Jost matrices) 



S\k) = S'\k) = S*{k) = S{-k), Im A; = , (4) 
lim S{k) = I, (5) 

Fl{k) = F+{-k) = F4k) , lm/c = 0, (6) 

/ is the unity matrix. As usual the physical formulation of the problem 
requires, in addition, the asymptotic condition 

hrn F+ik) = I (7) 
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and the condition 

detF+(mj) = 0, Kj > ,j = l,2, . . . ,m , m < +00 , (8) 
to be satisfied (m is the number of bound states). 

II Boundary value problem in the 
matrix case 

In the scalar (non-matrix) case of uncoupled partial waves the solution 
of the problem formulated above can be given explicitly []T7j]: 

F.ik) = n,(fc) exp - / » dk'\ , (9) 

n±(fc) = n||^, (10) 

U±{k) = l, m = 0. (11) 
However, in general case the matrices given by the Eqs.(0)-([TTD do 



not satisfy the boundary value condition (|ID • In fact, using the well- 
known equation 

1 1 ^ 

± m5{k' — k) 



k' -kTiO k' -k 
and substituting (0) in ([T|), one obtains 



^\g{k)+h{k) ^ ^g{k) ^-lg{k)+h{k) 



with 



-too ( h,l\ 

g{k) = HS-\k)Iil{k)) , h{k) = P.V. J 1^ dk\ 

— (X) 

P.V. stays for the principal value of the integral. So if we define the 
matrices A and B by 

A = g{k), B = -^g{k)+h{k) 



then the rule 

is to be fulfilled. The Eq. (|T2|) is vahd only if 

S{ki) S{k2) = S{k2) S{ki) , kij^k2. (13) 

The condition ( [131) means, from the physical point of view, that the 
potential in Schroedinger equation is central. If the interaction is of 
tensor character, then the condition (p^ ) and, consequently, the Eq. 
(|12|), are not fulfilled. This means that (0) can not be used for solving 
the boundary value problem of Sec. |. 

To construct the Jost matrix in this case let us begin with the first 
stage — preliminary factorization of the iS-matrix — "prefactoriza- 
tion" . 



Ill The S-matrix prefactorization 

Now the problem of finding of Jost matrix in the coupled channels 
case is reduced to the boundary value problem (|I]) with conditions 
(0)~(§)- 'S'-matrix can be diagonahzed by the use of an orthogonal 
transformation 

S{k) = U{k')[^ ^m,ik) j U-\k') , ImA: = 0, (14) 

where 5j{k) , j = l, 2 are the real-valued phase shifts, 

S^{-k) = - 6j{k) , ^(±oo) = 0, j = l,2. 
The matrix U (/c^) as usually can be written in the form: 



. 2\ _ ( cose(/c ) sm£{k ) 
\—sme{P) cos£(/c^) 

where £(/c^) is the real- valued mixing angle. We shall approximate 
the tangent of mixing angle e(/c^) by a rational function: 

tan.(.^) = gg. (15) 
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(16) 



Here and Qiv(^^) are polynomials in k of the degree 2M and 

2N, respectively. Without loss of generality one can suppose that the 
polynomials have no zeros simultaneously. The diagonal part of the 
^-matrix (|T^) can be factorized: 

e^^^^W \ / f,4k) \ / f,-_l{k) 
e^*'^^^^) ) [ f2-{k) )[ 

where fj±{k),j = l,2 are the solutions 

of scalar boundary value problems given by the Eqs.(|9|)- (pr[| ). Let us 
note, that if the number of bound states in the first channel (for j = 1) 
is equal to m\ then one must change m for m\ in (|9|)-([TT|). Similarly 
one has to change m for 777-2 for j = 2; nii + 7772 = 777 . Taking into 
account (|T5|) , (|T^) the problem ([I|) can be written in the form: 

F^^\k) = S-\k) F^^\k) (17) 

with 

So the preliminary factorization of iS-matrix is realized. 



IV The analytical properties of 
the matrix F^\k) 

The matrix F^\k) has a number of properties which are charac- 
teristic for the Jost matrix. As it follows from Eq. (|T^) F^\k) is 
analytical in the upper half-plane and the following equations take 
place: 

Fj°)*(/c) = Fi^\-k) = F^^\k) , ImA: = 0. (19) 

Nevertheless, F^\k) is not a solution to the boundary value problem 
of Sec. because: 1) the condition (|7|) is not satisfied; 2) det F^\k) 
has nonphysical extra zeros originated by the roots of the equations: 

QN{k') + iPM{k^) = 0, (20) 
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(21) 



Using ( P^ ) one can write 



so that the condition (g) is not satisfied. 

Let us emphasize one particularly important feature of the roots 
of the Eqs. (PI3|)-([?T|) . The polynomials are even functions of k and 
are not equal to zeros simultaneously (see Sec. |IT|) , so that the Eqs. 
( PD|) , (^) have no roots if Im /c = , or if Re /c = 0. Let us also no- 
tice that for real k the Eqs.(|2D|) , ( PT|) are complex conjugate. Now 
we can denote the roots of the Eqs. (|2I]|) and (|2T|) by ±iAj,±zA*, 
j = 1,2, . . . , L, L = max{M , A^} . For definiteness let us suppose 
Im(iAj) > , j = 1, 2, . . . , L . It is easy to see that the positions of all 
roots are symmetrical relatively to real and imaginary axes. We shall 
refer to these roots as to nonphysical zeros of the F^\k) -matrix de- 
terminant. The symmetry described in this Sec. is of great importance 
for the solution of the boundary value problem. 

Let us consider first the case of simple zeros. In order to obtain the 
Jost matrix, one needs only to remove the nonphysical zeros from the 
det F^\k) , but in such a way that the analytical properties of the 
matrix Fi'\k) remain unchanged and the condition ([7|) is fulfilled. 
Let us multiply the boundary value condition ([T7|) from the right by 
a rational real matrix W{k'^) 



If the matrix W{k'^) is such, that conditions on F\_ [k) mentioned 
above are fulfilled, then the Jost matrix F+{k) is of the form 



We will refer to the matrix W{k'^) as the correcting matrix. The 
method of obtaining the explicit form of W{k'^) is based on the ideas 



Fi^\k) Wi^) = S-\k) F^^\k) W{e) . 




(22) 



of H, M- 
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V The construction of the correcting matrix W{k'^) 

Let us make use of the symmetry property of the nonphysical zeros 
(Sec. |rV]) and let us fix four of them: ±iXi,±iXl . We shall now 
remove these nonphysical zeros from det F^^ (k) . To do this let us 
consider first the matrix 

where Pi is a projection matrix 

Pi = Pi. (24) 
Using the fact that (^ implies 

det Pi = , Sp Pi = 1 , 
it is easy to check that 

det^i(^-^) = (,^,.;,^,;.) - (25) 

Now we can obtain a new matrix: 

Fi'\k) = Fi^\k)Wi{e). 

As can be seen from Eq.(|2^) the det F^\k) has no zeros at the points 
k = ±^Ai, zbiAi- 

Now we need to verify the analyticity of F^\k) in the upper half- 
plane and the properties analogous to (|T^. This second condition will 
be satisfied if 

P* = / - Pi . (26) 

The explicit form of Pi must be chosen in such a way that the matrix 
F\^\k) will have no poles at the points k = iXi, iXl . Taking into 
account the Eqs. (^) , (^ , we can define 

Pi = ^yu (27) 
Spli 
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with 

Y - ( /2+(^Ai) U 1 W /r+(^Ai) \ 
'~[ /i+(^Ai) ; I, z 1 A /2V(zAi) J • ^'^^ 

Now F^\k) has no zeros in the upper half-plane. To check this fact 
it is sufficient to show that the following equations are satisfied: 

Ff(zAi)Pi = 0, (29) 

Fi'\iXl) {I - P,) = . (30) 

This can be done using the explicit form of F^\k) and Pi ( (|TB|) and 
( P7|) , (^), respectively) and the fact that k = iXi, iXl are the roots of 
the Eqs. (^) , (p). 

After removing of 4(n — 1), n < L nonphysical zeros we obtain for 
Fi")(/c): 

Fi''\k)=Fi''-'\k)Wn{k'), (31) 

^'^^^'^ ^^"-{k- iK){k + iK) + " (/c - iK) {k + zA* ) ' ^^^^ 

where P^^ has properties analogous to the Eqs. (^), (^) and is of 
the form ^ 

^"^S^^"' ^^^^ 
Yn = W-},{-Xl) X ... X M/f H-A^) (34) 



X 



/2+(ai) W 1 U /r+(^Ai) 

/i+(zAi) i U 1 M /2V(zAi; 



x[iyi(-A2)x...xiy„_i(-A 



2m* 



It is easy to show that det F^\k) has no nonphysical zeros k = iXj, iX*, 

j = 1, 2, . . . , n, and the matrix F^^\k) is analytical in the upper half- 
plane. 
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When all nonphysical zeros are removed and L [L = max{M, A^^}) 
projection matrices are obtained, then taking into account the Eq. (|7|) 
one has for the correcting matrix: 

W{k'^) = Wi{k'^) X ••• X WLik'^)C, (35) 

N>M, 

? "ol- (36) 



c = 



1 

PM 



QN PM ^ N = M, 

-PM QN J 

where Pm^Qn are the highest degree coefficients in the polynomials 
Puik^), QN{k'^), respectively. The Eqs. (^ , (^-(^ give the so- 
lution of the boundary value problem of Sec. in the case of the 
approximation (|T5D for the mixing angle. 

The boundary value problem of Sec. | can be solved in the case 
of multiple nonphysical zeros, too. However, in this case, one cannot 
write simple explicit equations for the projection matrix, although 
the idea of the construction remains to be the same. Let us notice, 
however, that, as was mentioned above, Jost matrices are used usually 
as an intermediate step during the solution of physical problems, and 
it is very useful to have definite simple form of them. For this purpose 
it is more convenient to have simple zeros. This is always possible to 
be done remaining within the limits of the experimental error, keeping 
in mind that the polynomials Pm(^^), Qnik'^) give an approximation 
for tan£(/c^). 



VI The example 

As an example let us consider the elastic two-channel neutron-proton 
scattering with the phase shifts ^Si , ^Di in the effective radius ap- 
proximation. In this case the scattering phase shifts and mixing angle 
are of the form: 

" (/c-m)(/c-Z9^)' t^^^i^i-2;^2' 
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where k = iK, is the bound state point, so that at m=l the condition 
is to be taken into account. 
The ^-matrix is a rational function and has the form: 



S{k) = 



1 / 2x' k^ 
k^ + [-k^ 2x^ 



{k+iK){k+i(p) Q 
{k~iK){k—iip) 

1 



2x2 

k^ 2x2 



Let us perform the iS-matrix prefactorization in terms of matrices 
Ff (fc) (see O, (0)) 



FS\k) = 



2x2 k^ 
-k' 2x2 



k^iK Q 
k±iifi 

1 



In our case the polynomials Pm(^2), QN{k^) (see (|T5|)) are of the form: 

PM{k') = k" , Qa.(/c2) = 2x2 , 

with M = 1, = 0. The number of projection operators Pj , j = 
1,2,... is L , L = max{M, A^} (see Sec. [V]) . In our case L = 1, so 
to obtain the Jost matrix we need to construct only one projection 
matrix Pi. The nonphysical zeros of (leiFf\k) (the roots of Eqs. 
(gg) , (gT])) are A; = ±Ai, ±A]; , with A = x(l-0- 

Following (p7|), ( |2B| ) it is easy to obtain the projection matrix Pi in 
the form: 

1 



2r/ 



Kx^ + {x-i^f) r]-ixi^ + K) 

v = x^ + {x + ^){x- 1^) ■ 

Using the condition (|7|) one has for the Jost matrix: 

k — lK Q 
1 

Pi 



F4k) 



2x2 e 
-k' 2x2 



k+iip 





X 



^ 7 + (/ — Pi) iXTT^ 



F+Af 



F+Af 



Instead of the Pi matrix we can use the P matrix defined by 





1 



-1 
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Now 



2x' 



-2x 
k 



2 k—iK 
k+iip 
2 k —iK 

-itp 



P 



/c2 + A? 



2* 



(37) 



This simple case considered as an example was investigated in [jT^ and 
one can show that ( |57| ) is nothing else that the well known result of 
T^] for the Jost matrix, but in more compact form. Our projection 



matrix P is just P3 — one of the projection matrices of [jT^. However, 
the (S-matrix prefactorization and the use of the symmetry property 
for nonphysical zeros make our algebraic method more effective and 
allows one to apply it in more complicated cases. 



VII Conclusion 

We have developed an effective algebraic method of the exphcit fac- 
torization of the (S-matrix in terms of the Jost matrices in the case of 
tensor potential. The only approximation used is the rational form for 
the mixing angle tangent. One of the main advantages of the method 
is the fact that the number of projection matrices involved does not 
depend on the phase shifts approximations. The obtained Jost matri- 
ces have simple compact explicit form and can be used in a number 
of physical situations (see e.g. [g]). 

Our result, as a mathematical result, is of independent interest by 
its own rights, too. We have enlarged the class of unitary matrices 
for which the solution of the Riemann - Hilbert boundary value prob- 
lem can be obtained in a way avoiding systems of singular integral 
equations. 



References 

[1] R. G. Newton, Scattering Theory of Waves and Particles^ 2nd 
ed.. Springer Verlag, New York (1982). 

[2] R. G. Newton, Inverse Schrddinger Scattering in Three Dimen- 
sions, Springer Verlag, New York, (1989). 



13 



[3] K. Chadan and P. C. Sabatier, Inverse Problem in Quantum Scat- 
tering Theory, 2nd ed., Springer Verlag, New York (1989). 

[4] H. V. von Geramb and H. Kohlhoff, Lecture Notes in Physics 
427, 285 (1994). 

[5] V. E. Troitsky, Lecture Notes in Physics 427, 50 (1994). 

[6] S.P.Novikov, S.V.Manakov, L.P.Pitaevskii and V.E.Zakharov, 
Theory of Solitons, Consultants Bureau, New York (1984). 

[7] R.K.Dodd, J.C.Eilbeck, J.D. Gibbon, H.C.Morris, Solitons and 
Nonlinear Wave Equations, Academic Press, New York (1982). 

[8] T. Fulton and R.G. Newton, Nuovo Cim. 3, 677 (1956). 

[9] V.E. Troitsky, Theor.Math.Fiz. 37, 243 (1978). [English transla- 
tion: Theor.Math.Phys. 37 (1978)]. 

[10] V.M. Muzafarov, Yad.Fiz. 27, 1686 (1978). [English translation: 
Sov. J. Nucl. Phys. 27 (1978)]. 

[11] R.G. Newton, J.Math.Phys. 21, 493 (1980). 

[12] V.M. Muzafarov, Theor.Math.Fiz. 58, 184 (1984). [English trans- 
lation: Theor.Math.Phys. 58 (1984)]. 

[13] R.G. Newton, J.Math.Phys. 31 10 (1990). 

[14] M.Coz, Lecture Notes in Physics 427, 353 (1994). 

[15] N. P. Vekua, Systems of Singular Integral Equations, P. Noord- 
hoff, Groningen (1967). 

[16] O.I.Kisaev, A.F.Krutov, D.I.Muravyev and V.E. Troitsky, Lecture 
Notes in Physics 427 163 (1994). 

[17] F. D. Gakhov, Boundary Value Problems, Pergamon Press, Ox- 
ford (1966). 

[18] G.N.Chebotarev , Uch. Zapiski Kazan. Univ., Mat. 116 31 (1956) 
(in Russian). 

14 



[19] R.G. Newton and T. Fulton, Phys.Rev.107, 1103 (1957). 



15 



